arXiv:1506.03856v2 [cond-mat.other] 29 Oct 2015 


Generalized Efimov effect in one dimension 

Sergej Moroz, 1,2, 0 Jose P- D’Incao, l ' 3 '@ and Dmitry S. Petrov 4 ] 3 ] 

1 Department of Physics, University of Colorado, Boulder, Colorado 80309, USA 
2 Center for Theory of Quantum Matter, University of Colorado, Boulder, Colorado 80309, USA 
2 JILA, University of Colorado and NIST, Boulder, Colorado 80309-0440, USA 
4 LPTMS, CNRS, Univ. Paris-Sud, Universite Paris-Saclay, 91405 Orsay, France 

We study a one-dimensional quantum problem of two particles interacting with a third one via a scale- 
invariant subcritically attractive inverse square potential, which can be realized, for example, in a mixture of 
dipoles and charges confined to one dimension. We find that above a critical mass ratio, this version of the 
Calogero problem exhibits the generalized Efimov effect, the emergence of discrete scale invariance manifested 
by a geometric series of three-body bound states with an accumulation point at zero energy. 


In quantum mechanics three identical bosons in three di¬ 
mensions interacting resonantly via a short-range two-body 
potential have an infinite tower of bound states, whose en¬ 
ergy spectrum forms a geometric series near the accumula¬ 
tion point at zero energy. This was discovered theoretically 
by Vitaly Efimov in 1970 |Tj] and is known today as the Efi¬ 
mov effect. This effect is a beautiful example of few-body 
universality since it is independent of the detailed form of 
the interaction potential provided it is tuned to the resonance 
(i.e., whenever a zero-energy s-wave two-body bound state 
if formed). The Efimov effect has been extended to systems 
of distinguishable particles |@-[6], liberated from three dimen¬ 
sions (7)] and found in other system s Qjj lgil. During the last 
decade a number of ex peri ments 110l4l4ll with cold atoms near 
Feshbach resonances 115(1 verified various universal aspects 
related to Efimov physics— the Efimov 4 He trimer has also 
been recently observed in JB— and demonstrated the exper¬ 
imental capability to explore fundamental aspects of few-body 
systems in exotic regimes. 

From a more general perspective, the most startling fea¬ 
ture of the Efimov effect is discrete scale invariance of the 
three-body problem, manifested in both bound and scattering 
three-body observables, that originates from continuous scale 
invariance of the two-body interaction. It thus appears nat¬ 
ural to us to generalize the Efimov effect to systems whose 
two-body interaction is not necessarily short-range and define 
it as the emergence of discrete scaling symmetry in a three- 
body problem if the particles attract each other via a two-body 
scale invariant potential 0. 

Motivated by this broader perspective on the Efimov effect, 
we study a three-body problem with a two-body long-range 
attractive potential of the form 


V(r) = 


a 

2 pr 2 


( 1 ) 


with /r being the reduced mass and a the dimensionless cou¬ 
pling constant. The potential 0 is scale invariant and, at 
zero energy or for sufficiently small r, where the energy term 
can be neglected, in one dimension the two independent so¬ 
lutions of the two-body Schrodinger equation are the powers 
r t/2±\A/ 4 -o^ However, the (inevitable) breakdown of the 
1 /r 2 law at small distances introduces a length scale b, made 


explicit by writing the linear combination of the two asymp¬ 
totic solutions in the form 


/ r\ l/ 2 +\/l/4—a / f \ 1 / 2 — \Jl /4 —a 

\b) ~ \b) 


( 2 ) 


For the further discussion it is crucial whether a is larger 
or smaller than 1/4 [4, TH B- The case a > 1/4 corre¬ 
sponds to the fall of a particle to the center and the discrete 
scaling is manifest already in the two-body problem. Here, the 
exponents 1/2 ± \/l/4 — a are complex conjugate, the two 
terms in Eq. ({2} should be treated on equal footing, and b be¬ 
comes an essential parameter, which can never be neglected. 
By contrast, the case a < 1/4 has two scale-invariant lim¬ 
its b = 0 and 6 _1 = 0 where, respectively, only the plus- 
branch r 1/2+ vT/ 4 -“ or only the minus-branch r V2-v' 1 / 4 -« 
survives in Eq. ( 0 . In practice, these two limits require, re¬ 
spectively, |6| -C £ or |6| f, where £ is a typical length- 
scale in the problem such as the system size, de Broglie wave 
length, etc. For instance, in Eq. 0 the minus-branch solution 
can be neglected if (fe/^)' 31-4 " <C 1. Thus, the plus-branch 
scaling is realized “automatically” by increasing the typical 
size of the system, whereas the minus-branch requires a fine 
tuning of the short-range part of the potential |20 h 22II . Phys¬ 
ically, this fine tuning signals the appearance of an additional 
two-body bound state emerging from the zero-energy thresh¬ 
old which can be realized using, for example, the Feshbach 
resonance technique B15ll . 


As far as the three-body problem with the two-body inter¬ 
action 0 is concerned, Calogero solved it in one dimension 
analytically for three identical particles 0 and found contin¬ 
uous scale invariance for all a < 1/4, which implies the ab¬ 
sence of the Efimov effect IB . In this paper we show that this 
conclusion does not hold in general for the modified Calogero 
problem - two identical spinless bosons or fermions interact¬ 
ing with a third particle via the potential 0 . In addition to 
the quantum statistics and the choice b = 0 or 6” 1 = 0 the 
modified problem is parametrized by the two continuous di¬ 
mensionless quantities: a and the mass ratio. Accordingly, 
we calculate the critical line separating the Efimov and scale- 
invariant regions and describe the nature of the three-body 
bound state spectrum in this parameter space. 
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The three-body Hamiltonian relevant for our problem reads 


where Ri and R 2 are the coordinates of two identical par¬ 
ticles of mass M and r is the coordinate of the third parti¬ 
cle of mass to. The potential V is given by Eq. CQ, where 
p = mM/ (M + to) and a denotes the interspecies dimen¬ 
sionless coupling. 

A convenient way to solve this problem is obtained using 
hyperspherical coordinates. First, we introduce the center-of- 
mass and mass-scaled Jacobi coordinates Rcm = \jnr + 
M{R ± + R 2 )]/(2M + to), x = yfpj2{2r - -Hi - R 2 ), 
y = ,/2ij,m(R2 ~ -Ri), where p = 2 mM/(m + 2 M) and 
Hm = M/2. It is then convenient to define polar (hyperspher¬ 
ical) coordinates x = 7Z cos 9, y = 7Z sin 9 with the mass- 
scaled hyperradius 7 Z = \/2'^ jj m l {r l — Rcm)' 2 - The inter¬ 
particle distances in the new coordinates become r — R\ = 
7Z sin (A + 9)/yJ2p, r — R 2 = Rsin(A — 6)/y/2p and 
R /2 — R t = 1Z sin O/y/M, where A = arctan \/l + 2 M/m. 
Accordingly, after separating the center-of-mass motion, the 
relative part of the Hamiltonian ([3]) is written as a two- 
dimensional radial problem 

H = -dn~ ^ M 2 e (4) 

with the hyperangular Schodinger operator 


Ml = ~dl - 


sin (A + 0) sin (A — 9) 


(5) 


9 = A 



FIG. 1. Hyperangular domain partitioning. 


solve the angular problem only in the domain 9 £ (0, 7t/ 2). 
Moreover, we will assume that the distinguishable particles 
are impenetrable. Physically, this is realized by regularizing 
the inverse square potential © with a short-range potential 
that has a strong repulsive core. Due to the interspecies im¬ 
penetrability, sectors I and II in Fig. |T] decouple and can be 
addressed separately. In sector I the hyperangular wave func¬ 
tion, i/j(9), should satisfy the following boundary conditions 
for 9 = 0 B 


ip = 0 fermions, 
ip' = 0 bosons 


(7) 


and for 9 —> A 


Two-body scale invariance leads to separability of the three- 
body problem in hyperspherical coordinates. The relative part 
of the three-body wave function ^>(TZ, 9) can thus be written 
in the factorized form ^(TZ, 9) = $(7Z)ip(9) and the problem 
separates into two tasks. First, one finds ip by diagonalizing 
the operator A4'/ 

Mgip = —s 2 ip. (6) 

Then, $(7 Z) is the solution corresponding to the Hamiltonian 
© with A ig substituted by — s 2 . This second task is trivially 
solved in terms of the Bessel functions J±i S , and the onset of 
the generalized Efimov effect coincides with the point s 2 = 0: 
for positive s 2 the system is Efimovian and for negative s 2 it is 
scale invariant. Thus, the problem of determining the critical 
mass ratio is equivalent to solving the hyperangular problem 
© and identifying the zero crossing of s 2 as a function of A. 
We will now discuss this procedure. 

The coincidence angles 9 = 0, tt (M — M coincidence) 
and 9 = ±A,7r ± A (M — m coincidences) partition the 
hyperangular circle into six regions (see Fig. 0}. Since two 
particles of mass M are identical, the wave function satisfies 
ip(6) = ip{—9) or ip(9) = —ip(—9), respectively, for bosons 
or fermions. In addition, the hyperspherical Hamiltonian is 
symmetric under 9 —> n — 9 and the wave function ip is either 
even or odd under this transformation. It is thus sufficient to 


ip ~ (A - 9) 1 / 2+VV4-C 
ip~ (A - 


plus-branch, 

minus-branch. 


( 8 ) 


The critical mass ratio is determined by solving Eq. © in 
sector I and is plotted in Fig. [2] (a) and (b) for bosons and 
fermions, respectively. We found that s 2 is an increasing func¬ 
tion of the mass ratio M/m for any choice of a and boundary 
conditions. In addition, we find no zero-energy (s = 0) solu¬ 
tion in sector II that satisfies the proper scale-invariant bound¬ 
ary conditions at the interspecies coincidence point. The wave 
function is thus zero in sector II, i.e., the probability to find 
the particle of mass to in between the two identical particles 
of mass M vanishes. 

It should be noted that the modified Calogero problem of 
the type © is exactly solvable and scale invariant for the plus- 
branch under the condition M/m = l/(l/2+ 1/4 — a) 

ESQ- The Efimov region corresponds to higher values of 
M/m and, since the problem is not solvable, we solve it nu¬ 
merically. We use the Numerov method [29]] on a logarithmic 
grid (see Ref. 12211 ). Nevertheless, we also find approximate 
analytic solutions for this problem in limiting cases discussed 
below. 

For the plus-branch the critical mass ratio diverges at a = 
1/8. In fact, for a > 1/8 both branchesgive rise to the 
Efimov effect for sufficiently large M/to| 30]. Indeed, for 
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FIG. 2. (Color online) (a) Critical mass ratio as a function of a for 
bosons: upper blue (lower orange) line is the Numerov numerics for 
the plus (minus) branch, dashed lines are analytic asymptotes near 
a = 0 and a = 1/8 (see text). The shaded regions (i), (ii), and (iii), 
denote the regimes in which the Efimov effect does not occur, occurs 
for the minus-branch solution only, and occurs for both minus- and 
plus-branch solution, (b) same as (a) but for fermions. 


M/m —> oo the angle A = 7r/2 and the hyperangular po¬ 
tential in Eq. (|3 reduces to —2ct/cos 2 9. One can see that 
the hyperangular problem becomes Efimovian for 2 a > 1/4 
independent of the quantum statistics of the heavy particles 
and the branch choice. This means that the spectrum of A i 2 g 
is unbound from below with deep bound states localized close 
to 6 = 7t/ 2. As a result, a finite 7r/2 — A is necessary to 
renormalize this potential and bring the ground state energy 
-s 2 to zero. Quantitatively, for the plus-branch solution in 
the vicinity of a — 1/8 we obtain II22[1 


— A w jYe 


±■ 5 -- 


(9) 


where the upper (lower) sign corresponds to the case of bosons 
(fermions) and jY = 16exp[—2 — 2\/2 — if/_ 3+v / 2 )/ 2 ] ~ 
11.887 with H n being the harmonic number. 

For the minus-branch the spectrum is Efimovian for any 
0 < a < 1/4 for M/m 1 1221 ■ The less stringent condi¬ 
tion for the Efimov effect in this case can be explained by the 
fact that the minus-branch two-body interaction nearly binds 


two particles and is, in this sense, more attractive than the 
plus-branch interaction with the same a. In fact, the hyper¬ 
angular problem can be solved analytically close to the non¬ 
interacting point a = 0. In Ref. (220 we show that for the 
bosonic case 


and for fermions 




( 10 ) 


( 11 ) 


The asymptotes (ITOl) . and (flTT) are plotted in Figs. [2] as 
dashed lines. 

The identified critical mass ratio is calculated using the 
wave function %l>(6) without nodes inside sector I. As one in¬ 
creases the mass ratio, wave functions with increasing number 
of nodes will give rise to additional towers of Efimov states. 

Now we describe the qualitative nature of the three-body 
bound state spectrum.The interaction in Eq. ([TJ must be reg¬ 
ularized at short distances, see 12211 . As the short-range po¬ 
tential depth Dq changes one can tune between a pure plus- 
branch (6 = 0) and minus-branch (6 _1 = 0) solutions. The 
nature of the three-body spectrum will depend on which re¬ 
gion in Fig. [2]the system falls into. There are three different 
regimes: 


• In the region (i), below the orange curve, there is no 
Efimov effect for any value of b. 


• In the region (ii), between the orange and blue critical 
curves, the spectrum behaves similar to the original Efi¬ 
mov problem [[TJ] . By starting from the plus branch so¬ 
lution with 6 = 0 and increasing the depth Dq, three- 
body bound states emerge one-by-one from the three- 
body continuum as one approaches 6 _1 =0. At the 
critical point 6” 1 = 0, where a zero-energy two-body 
bound state pops up, an infinite tower of Efimov states 
is formed with the Efimov parameter s_ (encoding the 
geometric factor e 27T/ ' s - f or the energy spectrum) which 
depends on both M/m and a. As one further increases 
the depth Dq, the trimers disappear one-by-one into the 
particle-dimer continuum. 


• In the region (iii) the spectrum resembles the one ap¬ 
pearing in the three-dimensional Efimov problem of 
particles with unequal scattering lengths 12l l3ll 1321. 


Now both the plus- and minus-branches support Efi¬ 
mov states characterized by the Efimov parameter s+ 
and S-, respectively, where 0 < s+ < s_. The 
energy spectrum contains an infinite number of three- 
body bound states close to the zero-energy threshold 
for any value of 6. The interpolation from the plus- to 
the minus-branch can be understood as follows: Near 
6 = 0 the energy spectrum close to the zero-energy 
threshold is controlled by the s + parameter. As one 
approaches 6 _1 = 0 the virtual dimer state of size of 
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^3 - ^3 -♦ 

R 2 Ri r 


FIG. 3. Two identical dipoles and a charge confined in one dimen¬ 
sion. 


order |6| is formed. As a result, the trimers with ener¬ 
gies below (above) ~ — 1//j b 2 follow the geometric 
scaling with the Efimov parameter s_ (s+). At reso¬ 
nance b ~ 1 = 0 the geometric spectrum with s_ scaling 
is obtained. When the depth Do is increased further, a 
two-body bound state is formed and three-body states 
with energies below (above) ed will also follow the ge¬ 
ometric scaling with the Efimov parameter s_ (s+) to 
the point where, when away enough from b~ 1 =0, the 
energy spectrum is again completely controlled by the 
s + parameter. 


Can the Efimov effect found in this paper be discovered 
experimentally? A promising candidate might be a mixture 
of dipoles and charges that are confined to one dimension. 
Indeed, the dipole-charge interaction in three dimensions is 
given by the scale-invariant anisotropic potential V(r) ~ 
cos 4>/r 2 1 where <j> is the angle between the direction of the 
dipole moment and the dipole-charge separation vector. Con¬ 
sider now a system of two identical dipoles of mass M and 
dipole moment P and a particle of mass m and charge q con¬ 
fined to a one-dimensional line. Let us also regard the dipoles 
as dumbbells with fixed dipole moments as illustrated in Fig. 
[3] This three-body problem is governed by the Hamiltonian 


H = -2M & 1 


-K e qP 






l(r-R r) 2 (r-i? 2 ) 2 J (R 1 -R 2 ) 3 ’ 


2K P P 2 


( 12 ) 


dipole-charge 


dipole-dipole 


where the Coulomb constant K e = l/(4tve). If we neglect the 
dipole-dipole term, this Hamiltonian maps on Eq. (0 with 


a = 2 fj.K e Pq (13) 

and thus gives rise to the Efimov effect provided a < 1/4 and 
the mass ratio is above the critical value. The presence of the 
dipole-dipole term introduces a length scale Idd ~ MK e P 2 ~ 
and bound states of this size. In the Efimov regime the 
length Idd provides the high-energy cutoff for the energy spec¬ 
trum. This cutoff should not affect the (geometric) energy 
spectrum close to the zero-energy threshold. However, since 
our problem is one-dimensional, the dipole-dipole interaction 
effectively fermionizes the dipoles since their wave function 
is suppressed at R\ — R 2 ~ Idd ■ Thus, the critical mass ratio 
for the Efimov effect should be read off of Fig. El (b) rather 
than (a) even if the dipoles are identical bosons. 


As an example, consider two polar molecules interacting 
with an electron. From Eq. (IT3l > the dipole moment of 
the polar molecule should satisfy P < P cr = eao/8 ~ 
0.318 Debye, where e is the charge of the electron and a o 
is the Bohr radius. Such a system has a large mass ratio, 
and, therefore, provided it falls into the region (iii) in Fig. 
|2](b), displays Efimov states (without fine-tuning to the mi¬ 
nus branch) which can be detected spectroscopically. For 
a typical mass ratio M/m = 10 5 we find s+ > 1 for 
P > 0.281 Debye. Moreover, Efimov states could also be 
observed if tuning of the dipole moment P is possible. In that 
case, near an electron-dipole resonance, dipolar losses should 
be enhanced every time a new Efimov state is formed. 


A three-dimensional version of the problem may have bet¬ 
ter chances to be realized experimentally and should also 
have interesting quantum-chemistry implications. In that case 

Pr , 


1.63 Debye |33| 


We consider it as a promising 


project and leave it for future studies. 
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SUPPLEMENTAL MATERIAL: GENERALIZED EFIMOV 
EFFECT IN ONE DIMENSION 


Logarithmic grid and Langer substitution 


Square well regularization and resonance 


Consider two particles with the reduced mass //. interacting 
in one dimension via the regularized inverse square potenital 


V U = -2^' 

VO) = -Co = - 


r > r-o, 
7 


2a» o 


0 < r < r 0 , 


(SI) 


A numerical solution of the Schrodinger equation 

[dg - v{0)\ip = k 2 V> (S6) 

with a potential possessing an inverse square singularity re¬ 
quires a grid which is dense near the singularity allowing to 
resolve the behavior of the wave function close to it. This can 
be accomplished by the Langer substitution 

6(x) = exp(x), tp{9) = exp (|)y(x) (S7) 


where ro > 0, 7 > 0. We will also assume V(—r) = V(r), 
which implies that a wave function should be either even or 
odd under the reflection r —> —r. It should be noted that 
the regularized potential necessarily supports a parity-even 
ground state with energy eos ~ — l/(/Li7g). This is in agree¬ 
ment with a general theorem that any attractive well supports 
at least one bound state in one dimension. Here, however, are 
interested in the zero energy solution. 

For r > ro a general solution of the relative Schrodinger 
equation at zero energy is given by 


which transforms Eq. < lS 6 b into 

dly{x) + g{x)y(x) = 0 (S 8 ) 

with 

g(x) = — exp(2x) |k 2 + i>[exp(x)] j — -. (S9) 

This equation is now solved by the Numerov method on the 
equidistant .x-grid with a lower cutoff ;j: mm < 0 . 


r ^ 1 / 2 +^ 1 / 4 -a ^ ( r , 1 / 2 -^ 174-0 


,, , , ( r 

ip{r) = A + [ — \ + 

Vr 0 / 


A-(±) 

\r 0 




(S2) 


On the other hand, for r < ro the solution at zero energy is 


ip{r) = C sin (y/7^) 
V’O") = Ceos {Vly) 


reflection odd. 


reflection even. 


(S3) 


Analytic calculation near a = 1/8 

Slightly above a = 1/8 for the plus-branch the critical 
mass ratio is large (A —> 7 r/ 2 ) and the angular problem can be 
attacked analytically. Here we sketch this calculation. First, 
it is convenient to redefine the angle to be (f> = A — 6. The 
angular Schrodinger equation ([ 6 } to be solved in sector I takes 
the form 


By matching the logarithmic derivative at r = ro one finds for 
the reflection odd case 



a 

sin 2 (</ + 2e) 


V> = o, 


(S10) 


Y=^± = V 1 / 4 - a - 1/2 + y/l cot ^7 rs;zn 
4 - VV 4 -a + i/ 2 -VrcotVr' 

For the reflection even case one must replace cot y/j —> 
— tan ^7 in Eq. (IS4b . For a generic value of 7 , the plus- 
branch is dominant for r > ro since X ^ 0. Note however 
that one can fine tune the regularized potential such that the 
wave function follows the minus-branch for r > ro. This 
happens for X = 0 resulting in the resonance condition 


Vy* — a — 1/2 + ypy cot -y /7 = 0 reflection odd, 

y/l/A— a — 1/2 — yfy tan ypy = 0 reflection even. 

(S5) 

By increasing the depth Dq of the short-range potential, the 
total potential supports more and more bound states. Inci¬ 
dentally, the resonance condition signals the appearance of an 
additional bound state emerging from the zero-energy thresh¬ 
old. 


where e = ir/2 - Ac 1 

In the regime cj)^> e this equation simplifies to 

(- 3 |- - r ^-)v , = o (Sll) 

V sin 


with a general solution 

^ =(1 — cos </>)i 


t pV 1 —8c*/2 


(cos (j>) 


L 2 Q^ 1 1/ 2°‘ /2 (cOS <j>) 


(S12) 


where P™(x) and Q™(x) are associated Legendre functions. 
The quantum statistics condition <0 imposed at 7 = A fixes 
the ratio L\/L 2 . 

In the regime </> <C 1 one can expand the denominators in 
Eq. (IS 1 01 1 and it transforms to 



a 

((f> + 2e) 2 


if) = 0 


(S13) 
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with a general solution 


Ip = ty/ x{2 + x) 

* + *) + S 2 < 3 ? S - 1 | /2 ( 1 + *) 


(S14) 


where x = (p/e. The ratio S1/S2 is now determined by apply¬ 
ing the plus-branch boundary condition ip ~ x 1 l 2 +\/ 1 l A - a a t 
x = 0 + [equivalent to Eq. ©]. 

In the intermediate regime x 1 and <p <C 1 the wave 
functions (IS 121 ) and (IS 1 4b must now be matched. In this re¬ 
gion they are given by 

ip = B+(/)^ 1+vl - Sa) + (S15) 

where B + /B_ = exp(i£) with a real phase £. The phases 
extracted from the solutions (IS 121 ) and (IS 14b can differ only 
by the angle 2nn, where n is an integer. The nodeless wave 
function ip is found for n = 1. The matching of the phases ( 
in the intermediate region with n = 1 gives rise to Eq. ©. 


In fact, the integrals on the right hand side of Eq. (IS 1 9b can be 
calculated for a constant ip, i.e., setting a = 0 in Eq. (IS 161) . 
Then, combining the result with Eq. (IS 1 7b we finally obtain 

- s 2 « 4a 2 - 2a{p- 1), (S21) 

from which Eq. (ITOl) follows immediately. 

For identical fermions the mass ratio diverges at a = 0. 
Hence one can follow the procedure described in the previous 
subsection with the wave function (IS 1 4b now satisfying the 
minus-branch boundary condition ip ~ x 1 / 2-1 / 1 / 4- " at x = 
0 + . As a result, one finds Eq. CD. 


Born-Oppenheimer approximation for minus-branch 

Here we study the minus branch in the regime M m, 
where the Born-Oppenheimer (BO) approximation can be 
used and the emergence of the Efimov effect can be intuitively 
understood. Within this approximation the three-body wave 
function is factorized 


Analytic calculation near a = 0 


1 '(R,r) = $(R)il>(r-,R), (S22) 


Here we show how the mass ratio for the minus-branch can 
be determined analytically close to the non-interacting point 
a = 0. In this case identical bosons and fermions must be 
treated differently. 

In the case of identical bosons the critical mass ratio van¬ 
ishes at a = 0. We can thus start from the analytical solution 
of the angular equation ([6]) at A = 7 t/4 


ip ex sin 1/2 -V 1 / 4 -“(6) + |) sin 1/2 -\ /1/4 -“(6» - (S16) 

with the positive eigenvalue 


- s 2 = 2(1 - 2^1/4 - a - 2a) = 4 a 2 + 0(a 3 ). (S17) 

Let us now take into account a small deviation of A from 7r/4 
perturbatively. In order to do this we redefine the angle 9 = 
4A0/-7T such that the configurational space is the fixed interval 
9 e (0, 7 t/ 4) independent of A. After this redefinition, Eq. ([6} 
becomes 



ap 2 

sin 2 [p(f +<?)] 


ap 2 

sin 2 [p(f-0)] 




-s 2 p 2 ip, 

(S18) 


where p = 4A/7T > 1. The linear shift of —s 2 (p) with respect 
to p — 1 <C 1 is now determined by the first-order perturbation 
formula 


d{~s 2 p 2 )/dp= {ip\v\\ip) / {ip\ip), (S19) 

where v\ is the first derivative of the potential in Eq. ( IS 181 ) 
with respect to p at p = 1 and ip is given by Eq. (IS 1 6b . We 
have explicitly 

vi(6) = — a sec 3 (29) 

x [ — 37r + 8 cos(2 9) + 7 r cos(40) + 160 sin(20)]. 

(S20) 


where the function ip(r ; R) corresponds to a bound state of the 
light particle in the combined potential of two heavy centers 
located at R\ = 0 and R 2 = —R. This wave function satisfies 


1 r 

2 TO - 


a 

I2 


a 


(•r + R) 2 - 


ip = eip. 


(S23) 


Hlight 


Note that the distance R is the only characteristic length scale 
in Eq. ( IS23I ). Therefore, if a bound state exists at some R, it 
also exists at any R, and its energy equals = —a/mR 2 , 
where u is a positive dimensionless number which depends 
on a and on the choice of the light-heavy boundary condition. 

Next, the Born-Oppenheimer approximation assumes that 
the light particle adiabatically follows displacements of the 
heavy ones for which acts as the effective interaction po¬ 
tential. Thus, the Schrodinger equation for the relative motion 
of the heavy particles reads 


d R + £GS (R) 


$ = E$. 


(S24) 


The energy spectrum E determines the energy spectrum of the 
three-body system. The spectrum of Eq. (IS24b is Efimovian, 
E n ~ where s ~ \JcrM / m. For simplicity, we 

neglected in Eq. (IS24I) the diagonal correction to the potential 
Vdia g (i?) = {' l P\dji' l P)/M which is subleading and does not 
modify the argument above. 

Now we demonstrate that Eq. (IS23I) actually supports at 
least one bound state for the minus-branch condition for any 
0 < a < 1/4. To this end, we consider the minus branch 
variational wave-function which equals 


ip K cx r 1 / 2 v /l / 4 a exp(— nr) 


(S25) 
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for r > 0 and vanishes for r < 0. We find that for any 0 < 
a < 1/4 and for a sufficiently small variational parameter k 

{ip K \Hi ight \ip K ) < 0. (S26) 

This can be seen by evaluating analytically the variational en¬ 


ergy (1S26I ) and Taylor expanding it around k = 0. We thus 
arrive at the conclusion that for the minus-branch boundary 
condition the BO approximation demonstrates that the en¬ 
ergy spectrum of our three-body system is Efimovian for any 
0 < a < 1/4 in the limit M m. 




